Effect of hybridization on the magnetic properties of correlated two-band metals 
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The magnetic properties of transition-like metals are discussed within the single site approxi- 
mation, which is a picture to take into account electron correlations. The metal is described by 
two hybridized bands one of which includes Coulomb correlation. The presented results indicate 
that ferromagnetism arises for adequate values of hybridization (V), correlation (U) and occupation 
number(n CT ). Some similarities with Dynamical Mean-Field Theory (DMFT) are indicated. 

PACS numbers: 71.10.-w, 71.10.Fd 



I. INTRODUCTION 

Recently [1 Q, the conventional view of the origin of 
ferromagnetism in metals has been under criticism. Tra- 
ditional mean-field calculations favor ferromagnetism but 
corrections tend to reduce the range of validity of that 
ground state 

Several recent works addressed the issue of ferromag- 
netism in metals, going beyond the Stoner model : Voll- 
hardt et al\^ furnish evidence of the stability of itinerant 
ferromagnetism in the one-band Hubbard model(HM)at 
electronic densities not too close to half-filling and large 
enough U. 

Nolting and BorgielQ use, for the one-band Hubbard 
model(HM), a spectral density approach(SDA), a two- 
pole ansatz that is equivalent to the Roth method^], 
containing some free parameters to be determined by fit- 
ting some spectral moments. Ferro-and antiferromag- 
netic states are obtained. Ferromagnetic solution occurs 
only if the occupation number n exceeds a critical value 
and U exceeds a minimum value. The ferromagnetic or- 
der arises due to a shift of the f and J, quasiparticle sub- 
bands. 

JarrelQ uses a self-consistent quantum Monte-Carlo 
method for the d = oo one-band HM but did not find a 
ferromagnetic state for any filling, temperature or corre- 
lation. 

Herrmann and NoltingQ, also in the one-band HM, 
use the spectral density approach (SDA), and a modified 
alloy analogy(MAA) method. In the latter, the atomic 
levels of the fictitious alloy are the ones found with the 
SDA. In both, a ferromagnetic state develops, although 
in the MAA it happens in a rather restricted region of oc- 
cupation n. They point out that an imaginary part in the 
self-energy, present in the MAA method but not in the 
two-pole ansatz, causes important differences in the mag- 
netic order of the model. An improved version of the AA 
which incorporates the intersite magnetic correlations Q 
seems to indicate a ferromagnetic instability for large U 
and n ~ 0.6. 
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Schwieger and Nolting [l| consider a two-band HM, as 
we do, but for d — » oo and use both the SDA and the 
MAA methods and show that interband particle fluctua- 
tions cause a spin dependent band shift and a spin depen- 
dent broadeninig of the Hubbard bands. This shift sta- 
bilizes and the broadening destabilizes ferromagnetism. 
For U — 5 the critical temperature T c is plotted ver- 
sus hybridization; strong fluctuations between the bands 
tend to suppress ferromagnetic order. 

We use instead a single-site approximation (SSA)[l0|. 
In this formulation, only one site, the origin say, exhibits 
the full Coulomb interaction; the others are subject to a 
local field: the spin and energy dependent self-energy S CT . 
S CT is sclf-consistently determined by imposing the van- 
ishing of the scattering matrix, a condition formulated in 
the Coherent Potential Approximation(CPA) for alloys 
and impurity problems, to restore translation invariance. 
We will be working at temperature T = 0. 

The metal is described by two non-degenerate bands, 
a and b; a is a Hubbard-like narrow band with in-site 
interaction U, hybridized (V a b) with the second one, b, a 
broader uncorrelated band. 



II. THEORETICAL MODEL: THE 
SELF-CONSISTENT SSA 

The starting Hamiltonian in this work is 



+ Y,{V ab bta ja + V+afb ja ) , (1) 

where nf a — a^aia and a denotes spin. t^ b is the tunnel- 
ing amplitude between neighboring sites i and j , in each 
band and V a b the hybridization. In the single-site approx- 
imation [l(| one adopts the following effective Hamilto- 
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In this expression 
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The method then replaces a translationally invariant 
problem, as defined by (|Tj). by an impurity problem where 
only the origin incorporates the Coulomb interaction, 
the other sites being acted by the local (k-independent) 
field But the effective Hamiltonian @, describes an 
'impurity problem' in presence of Coulomb intra-atomic 
term and we have to resort to some approximation. 

We use the Green function method [HI; after some 
algebra one obtains for the a-band Green function 
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V a \(k) 
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is the recursion relation of the a band including hy- 
bridization and 

a t a (cos(k x a) + cos(k y a) + cos(k z a)) 
e k = ^ ' (5) 

is the respective recursion relation of the bare band. In 
this paper we use t a = 1 and A = 3, in arbitrary en- 
ergy units. All energy magnitudes are taken in units of 
t a , making them dimcnsionless. The bare a band width 
is then W = 2. For simplicity we use e b k = ae k + e s 
(homothetic bands), for the bare b band. e s is the shift 
between the center of the bands. From now on we take 
fcjCt — > ki,i — x,y,z and V a b — Vb a = V ^constant in- 
dependent of kj. The scattering T^-matrix in Eq.Q is 
given by 



T a {w,T, a ) = 



U < ng_ CT> - + S CT (w)(C/ - E°'(w;))F CT (u;) 

[1-(U- + S CT (w)F CT (w)] 

I 
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where 
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The self-energy E CT is complex and spin dependent-thus 
giving a spin dependent band shift. The vanishing of 
the T-matrix gives a self-consistent equation for the self- 
energy: 

£ CT H = U< ng_ ff > +S CT H(?7 - T, a (w))F a (w). (8) 



The procedure presented here has some similarity 
to the one used in the Dynamical Mean-Field Theory 
(DMFT)[l3[in the following sense: here the original lat- 
tice model with correlation in every site is replaced by an 
effective one, where the correlation is present only at the 
origin ('the impurity') but at the same dimension d = 3. 



III. NUMERICAL RESULTS AND 
CONCLUSIONS 



It is important to stress that Eq.© results from a con- 
figuration average caracteristic of a CPA approach and 
not from a dynamical constraint. 



The Green function GJJ° then becomes 



w - e% - E CT (w) 



while Gf a is 



with 
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G&(w) = G b k (w) + V'G b k (w)GZ(w)G b k (w) (10) 



G b k (w) 



(11) 



the Green function of the bare b band. 



In the numerical calculations one chooses the total 
number of electrons per site as being about n = 2.2, a lit- 
tle more than half-filling. We start with < (n.?) >= 0.52 
and < (wj) >= 0.45 in eq©, searching for a ferromag- 
netic solution in a less than half filled a-band. The dy- 
namics generated by V and U, then, redistributes the a 
and the b band electrons, mixing the up and down states 
eventually producing a magnetization. In ref.[l( how- 
ever, the density of the correlated band, rather than n, is 
kept fixed. For comparison, in figs |T]) and ([2]), we took 
n = 2.2, a = 1.5 and e s = 1.0 while in fig ©, a = 2.5, 
n = 2.0 c e s = 1.0. 

We found that, for a certain range of parameters U 
and V, a ferromagnetic state (FS) develops. We then 
investigate the effect of V in this state, for a given U in 
the weak (U < W), in the intermediate (U ~ W) and 
in the strong coupling (U > W) limits. In fig ([1} we 
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FIG. 1: (Color online) Magnetization m =< n? > - < nj > 
of the correlated a-band versus V in a regime f7/W < 1, 
namely, U = 1. 



_N 0.015- 
03 



CD 



0.1 02 0.3 0.4 0.5 0.6 

V 



In the intermediate limit, e.g., U — 2, U/W ~ 1, as 
shown in fig ([2]), ferromagnetism is also observed ; again 
increasing V the metal tends to a non-magnetic state. 

In fig© we exhibit the magnetization for a typical 
strong coupling situation, namely, (7 = 5. Wc notice 
that the decrease of the magnetization with V is now 
slower than for smaller U. 

We have shown that the present method is computa- 
tionally feasible, producing reliable and physically sen- 
sible results compatibles with the existing literature on 
this important subject, opening new insights concerning 
the possible U/W regimes. Further work, for T > 0, and 
extending the model to include a two sublattice system 
in order to describe possible antiferromagnetic states, is 
in progress. 
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FIG. 3: (Color online) Magnetization m =< n" > - < nj > 
of the correlated a-band versus V in a regime U/W > 1, 
namely [7 = 5. 



FIG. 2: (Color online) Magnetization m =< n° > — < nj > 
of the correlated a-band versus V in a regime U/W ~ 1, 
namely U = 2. 



display the magnetization m versus V for {7=1. For 
V > 0.3 the magnetization is already very small or zero. 
For a given U, the effect of the hybridization is such as 
to reduce or suppress the magnetization. 
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